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Abstract 

We study the feasibility of performing a ry disappearance long-baseline experiment using 
a neutrino beam coming from muon decays. The advantage of such a technique with respect 
to the production of neutrino beams from pions is that in a muon decay both muon and 
electron neutrinos are produced in the same quantity. In case of ry —> ry jS oscillations, the 
v e charged current (CC) events can be used as a control sample, to predict in situ the 
rates, thus reducing the systematics due to the knowledge of the neutrino flux, which is the 
main source of uncertainties for disappearance experiments. We consider as our neutrino 
target, the ICARUS detector in its final mass configuration of 4.8 kton. 


1 Introduction 

The most satisfactory explanation of the observed disappearance of atmospheric muon neu¬ 
trinos result, strongly confirmed by the latest observations of the Super-Kamiokande[Q and 
Soudanll |J collaborations, is v — v T (or — v s ) oscillations. 

The importance of this result requires however confirmation in a different environment, 
i.e. using neutrinos produced from accelerators. 

Several proposals have been made to test the oscillations using a long baseline experiment, 
in Europe I Si !§• in Japan and in the US ||. The first long-baseline experiment, 
K2K, will start data taking in January 1999. This experiment, based on a flux of muon 
neutrinos from the KEK laboratory to the Kamioka mine, is expected to test ry —+ vx 
disappearance reaching the region down to Am 2 = 3 x 10 _3 eU 2 at 90% C.L. after three 
years of running. 

Two of the proposed experiments for the CERN-LNGS long-baseline project (ICARUSQ 
and OPERA 0) will be able to observe the appearance of the r in ry —> v T oscillation, and 
will ultimately reach the sensitivity of A m 2 about 1.5 x 10~ 3 eU 2 at 90% C.L. ©• 

While the region for oscillations suggested by the SuperKamiokande result is quite precise 
in the determination of the mixing angle consistent with maximal mixing, the uncertainty in 
Am 2 indicated by the SuperKamiokande zenith angle distribution is quite large, and does 
not exclude the possibility of this parameter to be as small as 6 x 10 _4 eU 2 . For this reason 
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it is important to envisage an accelerator-based experiment able to cover the full parameter 
space. 

The probability for neutrino oscillations between two families is given by 

P(E V ) = sin 1 2 26» sin 2 (1.27Am 2 —) (1) 

E v 

where 9 is the mixing angle between the two neutrino flavours, and Am 2 (in eV 2 ) is the 
difference of the squares of the neutrino masses. In the region with full mixing and small 
probability, it can be approximated as 

P(E V ) « (1.27Am 2 L/Ev) 2 (2) 

Thus, using a baseline of 732 km and neutrino energies of a few GeV, in the region suggested 
by the atmospheric neutrinos the oscillation probability is of the order of a few per cent. 

The ultimate sensitivity at low Am 2 in the —> vx disappearance mode will be limited 
by the knowledge of the neutrino flux. In fact, since the oscillation probability shows a 1/E 2 
dependence, the oscillated events are all concentrated in the low-energy region, where it is 
very difficult to observe any shape distorsion due to the oscillation. So, the only visible effect 
is a deficit in the total number of neutrinos observed. The integrated oscillation probability 
in this case is 


JP(E)^(E)dE J^(E)/E 2 dE 

J^(E)dE J^(E)dE 


(3) 


showing a quadratic dependence on Am 2 . 

In long-baseline experiments performed using a pion beam, the problem of the flux uncer¬ 
tainty has been addressed by the use of a near detector, this solution having the disadvantage 
that the solid angle covered by the two detectors is not the same, thus requiring a non-trivial 
extrapolation from the near to the far fluxjn] [ p~2[ . 

In the last few years there has been growing interest in studying the technical feasibility 
and the impact on physics of a muon collider with high intensity beams p~3||. Muons decaying 
in this machine are a continuous source of both muon and electron-like neutrinos ]l3|] via the 
following relationjj: 


/r —s- e u e u^,. (4) 

Between production and detection, neutrinos oscillate with a certain probability. Let us 
consider three possible cases of mixing between two neutrino flavours: 

1. Pure —> v T oscillations: in this case there will be disappearance of neutrinos, 

with oscillation probability P, given by equation [I]. If we define as <!>(, the initial 
flux, the oscillated flux <J> _ will be: 

<I\,. = x (1 -P) (5) 

The v T produced by the oscillation will have a flux 

( K = K. x P (6) 

but will not be seen in the detector because of the kinematical suppression resulting 
from the r mass. The D e component of the beam remains the same, since in this case 
neutrino electrons do not oscillate. The v e beam, however, provides a way to predict 

1 1ri this paper we always refer to the decay of negative muons. The same considerations are obviously valid in 

the case of positive muons. 
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at the far detector the initial flux. In the muon rest frame, the distribution of muon 
neutrinos and electron antineutrinos is precisely predicted by the V-A theory: 

d 2 N 12r 2 

_|. a _ ( (3_ 2l ) + ( l_ 2x) p„cos«] (7) 

T“Tir [(1 ^ ) + <1 - l)p - cos< ' 1 (8) 

where x = 2 E v /m^, P ^ is the average polarization of the muon beam, and 6 is the 
angle between the momentum vector of the neutrino and the mean angle of the muon 
polarization. In a muon collider, the beam polarization can be carefully measured via 
the spectrum of the electrons from muon decay; thus the spectra of the two components 
of the neutrino beam can be known with good accuracy, and so the ratio of the fluxes 
reaching the far detector . 

2. Pure —> i/ s oscillations: the flux and spectrum are exactly the same as for 

the —> v T case, and also the v e beam is not affected. It is possible, however, to 

discriminate between —> v s and v tJ —> v T oscillations, since the ratio of neutral 

current to charged interactions is different in the two cases. 

3. Pure Vy —> v e oscillations: in this case, also electron antineutrinos would oscillate into 

with the same probability as the v tJ —> v e oscillation. If we call this probability P 
(given as above by equation ^j), the final flux will be composed of four neutrino types: 



= *i M x(i 

-P) 

(9) 


= fcxP 


(10) 


= f‘ xP 


(11) 


= x (i 

- p ) 

(12) 


The total detected flux <IV (J + + 4>„ e + 3>t/ e in this case would be the same, i.e. 

equal to the initial flux. The antineutrinos coming from the oscillation would produce 
leptons of opposite sign with respect to neutrinos from the beam, so a magnetic detector 
with charge discrimination capability would be able to study in detail this kind of 
processes, as explained in |l5| for a different choice of the oscillation parameters. If 
no magnetic detector is available, it is possible to use the difference between neutrinos 
and antineutrinos cross-sections to get some hint for this kind of process. Let us call 
R the ratio between neutrino and antineutrino cross sections: 


a(y , CC) 
o(V,CC) 


(13) 


The fluxes of z7 e and v^ are equal at production, but since the energy spectra are differ¬ 
ent, the number of neutrinos reaching the far detector is slightly different for the two 
flavours. We assume this effect to be negligible, and we set ~ = $ (deviations 

from this behaviour can anyway be computed for a more accurate estimation). Thus, 
the number of interactions for the different kind neutrino flavours would be (assuming 
no charge discrimination): 


N^E) = N Vlt + JV„- = A x (1 - P) + A/R x P (14) 

N e (E) = N Ve + Nr e = A x P + A x (1 - P)/R (15) 


where A is the product of the flux times the neutrino cross section. The ratio between 
the two flavours is still dependent on P: 


N ^ _ 1-P + P/R _ 2 -P 
~N~ e ~ P+ (1 - P)/R ~ 1 + P 


(16) 
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where we have assumed R = 2, independent of energy. The above equation allows 
to derive the oscillation probability, with no assumption in the initial neutrino flux. 
Without charge assignment the interpretation of an observed deficit in terms of zy,, 
disappearance or —► v e oscillations is not straightforward, and the ambiguity could 
be solved considering data from other experiments (as usually done in the interpretation 
of the SuperKamiokande results) or studying the modification in the ybj distribution 
due to the the antineutrinos produced in v,, —> v e oscillations. 


2 Event rates 

The muons are produced in the decay chain of hadrons behind an appropriate target and 
are subsequently captured, cooled, accelerated and stored into a ring where they are let to 
decay. 

The muon yield is determined by several parameters. If a pulsed proton source is assumed, 
the total number of muons is given by the total proton yield times the average pion yield 
per proton, times the muon yield per pion, so that the total number of muons accelerated is 
given by the following relation: 


Nfj, — fbunch X Np/bunch X Y n /p X Y^j^ X t 

We assume a bunch repetition frequency of 10 Hz, with 10 14 protons per bunch. Given a 
standard yearly running of 10 7 s, it corresponds to 10 22 protons per year. 

Based on results from muon collider studies[|T3), we assume that the yield is about 0.6 
pions per proton, half of which can produce useful muons, and 0.25 muons per pion survived 
cooling and acceleration, thus resulting in about 0.1 muons per proton. Overall, a total muon 
production above 10 21 muons per year can be achieved. Let us remark that the requirement 
on cooling should be less stringent than that for a high luminosity muon collider where 
overlap between two high density beams is required. We assume that the total muon flux 
can be improved by a factor 2 with respect to the figures quoted above. By running the 
experiment during 5 years, a total of 10 22 muons could be reached. 

We assume that the storage ring would be placed inside the proposed tunnel for the 
CERN-Gran Sasso neutrino facility (NGS)Jh| (see Figure 0)- An optimized geometrical 
configuration of the storage ring is that of a ’’cigare”, in which two long straight sections are 
closed by two small arcs in which the muons are strongly bent. The radius of curvature (in 
meters) of a charged particle in a uniform magnetic held is given by 

R=J^ 

0.3 B 

where p is the momentum normal to the magnetic held in GeV, and B is the value of the held 
in Tesla. The straight section would ht in the decay tunnel of the NGS, which is foreseen to 
have a diameter of about 3 meters and a length of about 1 km. Given this configuration, the 
losses in the arcs can be neglected, and about 50% of the muon decays are directed towards 
the target. 

Using conventional magnets for the dipoles of the arcs, and assuming a magnetic held 
of 2 Tesla, the radius of curvature would be about 10 meters, therefore some modihcation 
would be needed to the proposed target and hadron stop caverns to accomodate the arcs. If 
on the other hand superconducting magnets could be used (8 Tesla), the radius of curvature 
would be around 2.5 meters, and both arcs could ht in the caverns proposed in the NGS 
project. 

In this conhguration, half of the neutrinos would be emitted in the backward direction 
with respect to the Gran Sasso laboratory. Given the slope of the tunnel, they would reach 
the surface in a site located about 1.5 km from the muon accelerator, between the CERN- 
Meyrin and CERN-Prevessin laboratories. There, a detector could profit from the extremely 
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Muon energy (GeV) 

v IJt CC events 

v e CC events 

2.5 

560 

240 

3.5 

1540 

660 

4.5 

3280 

1410 

5.5 

6000 

2570 


Table 1: Charge current rates in the ICARUS detector for neutrinos produced by the decay of 
10 21 muons, at different muon beam energies. 


intense beam to provide a further cross-check on the neutrino flux. This location would act 
as a near detector, but in a much more favorable location. 

In the other direction, the neutrinos directed towards the Gran Sasso would travel the 
distance between the two laboratories (732 km), to be detected in one of the large apparatus 
installed there. 

We consider the ICARUS Q detector as a target, in its final mass configuration of 
4.8 ktons. We have there an excellent energy resolution for electrons, and a very good 
electron-muon separation capability. We purposely ignore the charge measurement, not 
considering the possibility of studying —> v e oscillations using the charge of the lepton as 
a discriminator. 

The total number of z/ ;j and v e charged current interactions in ICARUS as a function of 
the muon beam energy is shown in Figure |] and in Table [l], considering unpolarized muon 
beams. The energy spectra of the two types of neutrinos for different values of the muon 
energy are shown in Figure ||. We have assumed for simplicity that the neutrino cross sections 
scale linearly with energy even in the very low-energy part of the spectrum: 

a cc (v e ) = a cc (y^) = 0.67 x 10~ 38 x E(GeV) cm 2 

<r CG (i7 e ) = a cc (^) = 0.34 x 10" 38 x E(GeV) cm 2 

The relevant point is that for the energies considered here the cross section for electron and 
muon neutrinos are the same. 


3 Results 

We consider the results that could be achieved by the experimental setup described above, 
without near detector, in the case of pure —► v TtS oscillation. 

For the different values of muon beam energy, the oscillation probability as a function of 
the parameter Am 2 is shown in Figure || In the region of small probabilities, the dependence 
of this quantity on Am 2 is quadratic, as seen from equation |j. The probability has also 
a quadratic dependence on the neutrino energy, so the low-energy muons yield a higher 
oscillation probability with respect to the high-energy muons. In the A?n 2 region of interest, 
between 10~ 4 and 10” 3 eU 2 , the oscillation probability is between 1 and 10%, and a similar 
or better sensitivity is needed for observing the oscillation signal. 

If the energy shape of the oscillated events is not considered, and only a counting experi¬ 
ment is performed, as is the case for the small-Am 2 region, the statistical error is just given 
by the total number of muon neutrinos observed: 

A Np = _J_ 

N, v^' 

The total error on the oscillation probability has to account also for systematic effects. The 
main systematic effect comes from the knowledge of the flux. Without the near detector, this 
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is given in first approximation by the statistical error on the number of electron neutrinos 
detected in the far detector. Other sources of systematics are corrections due to the knowl¬ 
edge of the muon polarization, and from uncertainties on electron and muon identification 
efficiencies: 

A P AN U A N e „ 

~p~ ~ 

The contribution of the first two components of the total error, of pure statistical nature, are 
plotted in Figure |] as a function of the integrated flux. Even in the most favorable case, this 
error never gets below 2%, and it is assumed that the other sources of systematics, indicated 
by S in the above expression, are controlled to a better level, so they are neglected in the 
total error estimation. 

Being of only statistical nature, the total error decreases as the square root of the total 
flux. If a near detector is built, the statistical error on the number of v e detected becomes 
negligible, but then systematic effects due to the extrapolation from the near to the far 
detector become relevant. 

The error also shows a dependence on the inverse of the square of the beam energy, 
because of the combined effect of the increased neutrino cross section and of the larger 
number of neutrinos reaching the far detector as an effect of the Lorentz boost. Since, as 
seen before, the oscillation probability decreases with the square of the beam energy, at first 
approximation the sensitivity of the experiment is independent on the energy. 

From the curves shown, it is possible to see that for the low energy beam the oscillation 
probability is equal to the lcr error in the case of 10 22 muons for values of Am 2 around 
4 x lCT 4 eU 2 . 

The above discussion is based on the fact that at low oscillation probabilities the oscil¬ 
lation occurs in the low energy part of the spectrum, where statistics is limited. At higher 
values of Am 2 , however, more information can be extracted from the energy distribution 
of the detected neutrinos, in addition to the total flux, since in this case the characteristic 
oscillation pattern can be seen. As an example, the spectra of v e . without oscillations, 

with oscillations and their ratio are shown in Figures and ^ assuming value of Am 2 of 
2.2 x 10 _3 eU 2 and 1.0 x 10~ 2 eU 2 , and considering the smearing due to the energy resolution 
of the ICARUS detector. The oscillation patterns are clearly visible. 

All considerations done so far have been made in the assumption of maximal mixing. 
When smaller mixing angles are also considered, two-dimensional oscillation contours can 
be extracted. For the final result we consider the two cases of integrated fluxes: 10 21 and 
10 22 muons. The 90% C.L. contours are shown in Figure for a muon energy of 2.5 and 
5.5 GeV. It can be seen that, while in the “low” intensity scenario the covered area reaches 
values of Am 2 of about 10~ 3 eU 2 , for the high intensity case the sensitivity goes to about 
6 x 10~ 4 eU 2 , thus covering at 90% C.L. all the range indicated by the atmospheric neutrino 
results. 


4 Conclusions 

The SuperKamiokande result, allowing the possibility of very small A in 2 for neutrino os¬ 
cillations, pose strong challenges to the experiments willing to explore this region of the 
parameter space. The method presented here is based on disappearance of coming from 
a muon beam. Given the layout of the accelerator, a near position could be placed on the 
surface, opposite with respect to the far laboratory. Since a built-in control sample of z7 e 
coming from the same muons can be exploited for normalization, the experiment can how¬ 
ever be performed without requiring a near detector, and in this configuration it has the 
potentiality of reaching in a clean way the low-Am 2 region suggested by the atmospheric 
neutrino results. 
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Figure 1: Layout of the tunnel for the CERN-Gran Sasso beam 
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Figure 3: Energy spectra of (upper plot) and v e (lower plot) interacting via charged current 
processes for muon energies of 2.5, 3.5, 4.5 and 5.5 GeV (10 21 muons, 4.8 kton). 



















































%40f _I—I 

^30 3 - 

^ 20 Y - - 

10 ^ r — 

r\ Ej i i i l i i i i — 1 i i i i l i i i i l i i i i l i i i i l i i i i l i i i i l i i i i l i i i i 

0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 


v e events (ICARUS, 10 2 V) Ge V 



v (oscillated) events (ICARUS, 10 21 \i) GeV 


Figure 6: Energy spectra of detected u e [upper plot], non-oscillated (solid line) and oscillated 
(dots) [middle plot] and the ratio of the two [lower plot] for 10 21 muons of 5.5 GeV and 
Am 2 = 2.2 x 10 -3 eV 2 
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Figure 7: Energy spectra of detected i7 e [upper plot], non-oscillated (solid line) and oscillated 
v jx (dots) [middle plot] and the ratio of the two [lower plot] for 10 21 muons of 5.5 GeV and 
Am 2 = 1.0 x 10 -2 eV 2 
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Figure 8: 90% C.L contours for a disappearance experiment with a total flux of 10 21 and 10 22 
muons, compared to the SuperKamiokande 90% C.L. contour. Full line is the limit for 5.5 GeV 
muon beam and 10 22 muons, dashed bold line is for 2.5 GeV beam and 10 22 muons, dashed line 
is for 5.5 GeV beam and 10 21 muons, dot-dashed line is for 2.5 GeV beam and 10 21 muons. 
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